Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Initialising ...
Dehbi, A.*; Cheng, X.*; Liao, Y.*; 岡垣 百合亜; Pellegrini, M.*
Proceedings of 19th International Topical Meeting on Nuclear Reactor Thermal Hydraulics (NURETH-19) (Internet), 15 Pages, 2022/03
Nuclear degraded cores produce fission product aerosols that may reach the environment if not removed by natural processes and/or filtering equipment. The transport paths of aerosols usually include transits through stagnant water pools. It is therefore essential to develop computational tools to predict the aerosols retention by water pools. Currently, this is mostly done with 1-D lumped-parameter codes that are too simplistic to capture the physics. It is hence worthwhile to attempt the CFD approach, which has recently become reasonably mature to address bubble hydrodynamics in low momentum two-phase flows. In this first comparative exercise, we restrict the investigation to a hypothetical parallelepiped water pool (2
28 cm
) into which air is injected through a circular 4 mm ID orifice at low velocity of 0.2 m/s. We present predictions of the gas phase dynamics (void and velocity profiles) for both Euler-Euler and Interface Tracking (IT, Volume-of-Fluid (VOF)) methodologies. In addition, we compare bubble shape, volume and detachment frequency from various IT simulation codes (CFX, Fluent, Star-CCM+, OpenFOAM). Reasonable agreement is found between IT simulations near the injector, but discrepancies increase as one moves towards the free surface. The disagreement between the Euler-Euler and IT results is substantial throughout the domain. Future studies will consist of validation exercises against experimental data to highlight potential model deficiencies and point to ways of remedying them.
齋藤 慎平*; De Rosis, A.*; Fei, L.*; Luo, K. H.*; 海老原 健一; 金子 暁子*; 阿部 豊*
Physics of Fluids, 33(2), p.023307_1 - 023307_21, 2021/02
被引用回数：34 パーセンタイル：98.3(Mechanics)流れ場で沸騰が発生する現象は強制対流沸騰として知られている。今回、飽和条件の流れ場中の円柱上での沸騰システムを数値的に調査した。複雑な気液相変化現象に対処するために、擬ポテンシャル格子ボルツマン法(LBM)に基づく数値スキームを開発した。高いレイノルズ数の数値安定性を高めるため、衝突項を中心モーメント(CM)の空間で解いた。CMベースのLBMに適した力場スキームを採用することで簡潔でありながら堅牢なアルゴリズムとなっている。さらに、熱力学的一貫性を確保するために必要な追加項をCMの枠組みにおいて導出した。現在のスキームの有効性は、核形成,成長、および30-30000の間で変化するレイノルズ数の蒸気泡の離脱を含む一連の沸騰プロセスに対してテストされた。開発したCMベースのLBMは、初期気相などの人工的な入力なしに、核沸騰,遷移沸騰、および膜沸騰の沸騰様式を再現できる。結果からプール沸騰ではなく強制対流システムでも抜山曲線として知られる典型的な沸騰曲線が現れることが分かった。また、今回のシミュレーションは膜沸騰領域でも断続的な直接固液接触の実験的観察を支持することが分かった。
Wan, T.; 直江 崇; 粉川 広行; 二川 正敏; 大林 寛生; 佐々 敏信
Materials, 12(4), p.681_1 - 681_15, 2019/02
被引用回数：3 パーセンタイル：17.62(Chemistry, Physical)To perform basic R&D for future Accelerator-driven Systems (ADSs), Japan Proton Accelerator Research Complex (J-PARC) will construct an ADS target test facility. A Lead-Bismuth Eutectic (LBE) spallation target will be installed in the target test facility and bombarded by pulsed proton beams (250 kW, 400 MeV, 25 Hz, and 0.5 ms pulse duration). To realize the LBE spallation target, cavitation damage due to pressure changes in the liquid metal should be determined preliminarily because such damage is considered very critical from the viewpoint of target safety and lifetime. In this study, cavitation damage due to pressure waves caused by pulsed proton beam injection and turbulent liquid metal flow, were studied numerically from the viewpoint of single cavitation bubble dynamics. Specifically, the threshold of cavitation and effects of flow speed fluctuation on cavitation bubble dynamics in an orifice structure, were investigated in the present work. The results show that the LBE spallation target will not undergo cavitation damage under normal nominal operation conditions, mainly because of the long pulse duration of the pulsed proton beam and the low liquid metal flow velocity. Nevertheless, the possibility of occurrence of cavitation damage, in the orifice structure under certain extreme transient LBE flow conditions cannot be neglected.
小野 綾子; 鈴木 貴行*; 吉田 啓之
Proceedings of 26th International Conference on Nuclear Engineering (ICONE-26) (Internet), 6 Pages, 2018/07
液膜蒸発モデルに基づいた限界熱流束予測評価手法の開発において必要な伝熱面ごく近傍の気泡の挙動や気泡同士の接合現象について、詳細二相流解析コードTPFITによる発泡,気泡生長および接合をシミュレートすることで知見を得る。本発表では、高圧域での沸騰現象の理解と、気泡形状の影響の理解のために、圧力をパラメタとした解析結果に基づいて、気泡の挙動に圧力が与える影響について報告する。
二川 正敏; 直江 崇; Tsai, C.-C.*; 粉川 広行; 石倉 修一*; 池田 裕二郎; 祖山 均*; 伊達 秀文*
Journal of Nuclear Materials, 343(1-3), p.70 - 80, 2005/08
被引用回数：57 パーセンタイル：95.58(Materials Science, Multidisciplinary)核破砕中性子源水銀ターゲットには、大強度のパルス陽子ビームが入射される。このとき、水銀内には熱衝撃に起因する膨張波が生じ、その伝播過程でキャビテーションによる損傷がターゲット容器内壁に形成される。この損傷挙動を電磁力衝撃圧負荷試験装置により観察し、繰り返し周波数及び入力パワー依存性について調べた。さらにバブル崩壊により誘発される音響振動に着目した損傷ポテンシャルを提案するとともに、バブルダイナッミックスに基づく数値解析を行った。その結果、以下の結論を得た。(1)侵食挙動は微小領域塑性変形支配の潜伏期と質量減少を伴う定常期に大別できる。(2)定常期の質量減少を予測しうる実験式を導出した。(3)潜伏期間は、入力パワーの増加とともに縮小し、繰り返し周波数の増加により拡大する。
井田 真人
Journal of the Physical Society of Japan, 73(11), p.3026 - 3033, 2004/11
被引用回数：8 パーセンタイル：48.25(Physics, Multidisciplinary)音響的に相互作用する二気泡の「遷移周波数」について数値的に検証した。Ida[Phys. Lett. A 
(2002) 210]により、音場中の二気泡系に三つまでの遷移周波数が存在することが理論的に示された。これは気泡の脈動と入射波の位相差を
にする特徴周波数である。続編にあたる論文[M. Ida, Phys. Rev. E 
(2003) 056617]ではさらに、遷移周波数こそが二次Bjerknes力の符号を決定する要因であることが理論的に示唆された。本論文では表面張力項付きの圧縮性Navier-Stokes方程式を支配方程式に選んだ直接数値シミュレーションによってこれらの理論的結果を検証する。得られた数値計算結果はそれらの理論的予測をよく再現し、遷移周波数の存在を正当化する。
石倉 修一*; 粉川 広行; 二川 正敏; 神永 雅紀; 日野 竜太郎; 齊藤 正克*
日本原子力学会和文論文誌, 3(1), p.59 - 66, 2004/03
原研とKEKは共同で、大強度陽子加速器を用いた複合研究施設(J-PARC)の建設を進めている。その中核設備となるMW規模の核破砕中性子源に用いる液体水銀ターゲットが約1MWのパルス状陽子ビーム(パルス幅1ms)を受けた時、瞬時の熱膨張による圧力波が発生し、ターゲット容器に負荷すると同時に、キャビテーションの発生原因ともなる。キャビテーションが崩壊する時に誘起される水銀マイクロジェットが固体壁へ衝突する時のマイクロピット形成挙動を検討するために、水銀の状態方程式の非線形性と固体壁材料の降伏応力の歪速度硬化依存を考慮して、単一気泡を対象とした気泡動力学に基づき、水銀球が固体壁に衝突する解析を実施した。その結果、静的降伏応力が高い材料の場合、ピットの深さを抑制できるなどの、マイクロピットの形成と材料強度との関係を理解するうえで、重要な結果が得られた。
石倉 修一*; 粉川 広行; 二川 正敏; 菊地 賢司; 日野 竜太郎; 荒川 忠一
Journal of Nuclear Materials, 318, p.113 - 121, 2003/05
被引用回数：12 パーセンタイル：62.16(Materials Science, Multidisciplinary)出力1MW,パルス幅1
sの陽子ビームが入射した時に水銀ターゲット容器に発生する熱衝撃応力を解析した。初期圧52MPaが水銀中を伝播することにより61MPaの最大負圧を生じる。この負圧の発生によりキャビテーションの生成が予想される。そこで、キャビテーションの挙動を調べるために、単一バブルに対してバブル動力学から導かれる運動方程式のシミュレーションを行い、パラメータ解析を実施した結果、ターゲット容器ウィンドウ部の圧力が変化することによりバブルが1000倍以上に膨張することがわかった。その結果、波動伝播に影響を与えることになるが、バブル液体の状態方程式は多項式近似が可能である。また、静的・動的試験により液体の臨界圧力,音速,共振周波数を測定することにより、水銀中に内在するバブルの半径や体積率の判定が可能となることが示唆された。
功刀 資彰
日本機械学会論文集,B, 63(609), p.88 - 96, 1997/05
本研究は、VOF(Volume of Fluid)法の特徴である質量の完全な保存性を保ちつつ、その欠点である界面こう配輸送の不正確性を克服する新しい流体率輸送法を開発し、最終的には相変化を伴う固体・液体・気体の混在する多相熱流体の統一的直接解析手法の構築を目指している。本報告では、表面張力モデルの取り扱い及び新しく開発した「界面こう配を考慮した界面輸送法」を含む多相流の直接解析手法MARS(Multi-interface Advection and Reconstruction Solver)の概要を述べるとともに、解析例として、2次元気泡塔下部オリフィスから連続的に発生する気泡及び発生した気泡同士の合一挙動、3次元気泡塔内における2つの気泡の合一挙動及び2次元容器内における沸騰・凝縮挙動シミュレーションを示す。
